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A new method to constrain gravitational theories depending on the Ricci scalar is presented. It is 
based on the weak energy condition and yields limits on the parameters of a given theory through the 
current values of the derivatives of the scale factor of the Friedmann-Robertson- Walker geometry. A 
further constraint depending on the current value of the snap is also given. Actual constraints (and 
the corresponding error propagation analysis) are calculated for two examples, which show that the 
method is useful in limiting the possible f(R) theories. 



Introduction 

It follows from several observations Q that the uni- 
verse is currently expanding with positive acceleration. 
The many models that have been advanced to explain 
this situation can be classified in two classes. The first 
class contains those models that incorporate modifica- 
tions to the matter side of Einstein's equations. This 
matter (known as "dark energy"), can be described by 
an ideal fluid or by a scalar field with a convenient poten- 
tial Q, and it must violate the strong energy condition 
in order to accelerate the universe in General Relativity 
(GR) . The models in the second class have normal matter 
as a source but assume that gravitation is not described 
by GR at low curvatures. As examples of this latter class 
we can mention theories depending on the Ricci scalar 
S (the so-called f(R) theories), and gravity modi- 
fied from contributions of extra dimensions [5j . Although 
f(R) theories offer a chance to explain the acceleration 
of the universe, they arc not free of problems. For in- 
stance, the application of their metric formulation to the 
Friedmann-Robertson- Walker (FRW) geometry yields a 
fourth order nonlinear differential equation for the scale 
factor a(t) which in general cannot be analytically solved 
even for simple f(R). 

Since many f(R) give rise to accelerated expansion, 
another issue is how to reduce the theory space using ob- 
servations. Constrains have been obtained from cosmo- 
logical and astrophysical data, solar system tests, fifth 
force/BBN data, and by requiring that a given theory 
describes the correct sequence of decelerated-accelerated 
phases in the evolution of the universe . Most of the 
cosmological tests involve either some transformation of 
the theory under scrutiny to an equivalent form with one 
auxiliary scalar field and/or considerations in different 
frames (see for instance 0), or some assumptions regard- 
ing the dependence of the Hubble "constant" H with the 
redshift 81 . Here instead a new criterion based on model- 
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independent data shall be given, that helps in deciding, 
without solving the EOM or making frame transforma- 
tions or assumptions about H, whether a given f(R) the- 
ory is appropriate to describe the universe. The basic 
premise will be that the acceleration is due solely to a 
modified gravitational theory with normal matter as a 
source. The criterion will then be obtained by imposing 
the energy conditions on matter, yielding conditions on 
f(R) and its derivatives w.r.t R in terms of the current 
value of the derivatives of the scale factor. These con- 
ditions are to be satisfied if the theory given by f(R) 
is to describe the current state of the universe, and they 
bring forth limits on the parameters that enter the theory 
under consideration. 



Energy Conditions 

The energy conditions (EC) are inequalities satisfied by 
"normal" matter (see for instance |9j)- When specialized 
to a FRW universe, the (local) null, weak, strong and 
dominant EC are given by 

NEC p + p > 0, 
WEC p > and (p+p > 0), 
SEC<=> (p + 3p>0) and (p+p>0), 

DEC p > and (p±p > 0). 

The EC have proved to be useful in the context of cos- 
mological singularities ^(| an( A bounces Other ap- 
plications of the energy conditions to cosmology can be 
found in 0. 

Let us remind the reader that for a theory given by 
f(R), the EOM are 

f'R - 2/ + 3/" ( R + — ) + 3f"R 2 + T = 0, (1) 
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and we have assumed a flat universe. From Eqns.Q and 
© the energy density and the pressure of the fluid can be 
expressed in terms of the scale factor and its derivatives: 
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Before proceeding to build with these equations the in- 
equalities that define the energy conditions, let us remark 
that observations show that the current matter content 
of the universe (assumed here to be normal matter, as 
opposed to dark energy) is pressureless. In this case the 
EC reduce to the inequality po > plus the equation 
Po = 0, where the subindex means that the quantity 
is evaluated today. We shall express these conditions in 
terms of following kinematical parameters: the Hubble 
and deceleration parameters, the jerk, and the snap, re- 
spectively given by Q 
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While the current value of the first three parameters 
is H = 72 ± 8km/(secMpc) 0, go = -0.81 ± 0.14, 
jo = 2.16^0 75 16], no measurements of the snap have 
been reported yet. By writing po > in terms of the 
parameters we get 



3<ZoH 2 /o - y - 18#o 4 /o tio 



qo - 2) > 0. (5) 



This inequality gives a relation between the derivatives 
of a given f(R) and, as will be seen in the examples 
of the next section, it limits the possible values of the 
parameters of the theory. 

Notice that Eqn.Q involves the snap (through R). If 
we had a measurement of so, we could use the equation 
Po = to obtain another constraint on f[R). Since this 
is not the case, we shall express po — in such a way 
that it gives the possible current values of the snap for a 
given f(R): 
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Examples 

To see how Eqns.0 and |JBJ can be used to put con- 
straints on a given f(R), let us examine two examples. 
The first one is given by ^3] 



f(R) - aR- n , 



(7) 



with n £ A/", which can be taken as the n = 1 low- 
curvature limit of 



f(R) = R 



a 
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a model studied in §|. Substituting Eqn.Q) in Eqn.JSJ) 

for a > and n even we get 

-3qoH^nR -^R 2 ~18H^n{n + l)(jo-q -2) > 0. (8) 

Replacing in this equation the numerical values of the 
parameters and using Ro — 6Ho(qo — 1) we get an in- 
equality that must be satisfied by n: 



-17.64n 2 - 44.50n - 59.62 > 0. 



(9) 



Since this equation cannot satisfied by any real n, we 
conclude that n cannot be even for a > 0. The same 
analysis with odd n reverses the sign of the inequality 
©, so only odd values of n are allowed for a > 0. This 
result generalizes that obtained in Q for n = 1 |2l|. In 
the same way, we obtain that only even n are allowed 
for a < 0. These conclusions are valid even when the 
error coming from the kinematical parameters is taken 
into account (see Fig|lJ. If we knew the value of s , we 
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1: Plot of <f> (solid curve, see Eqn.Q) and the associated 
8(f> (dashed curve) in terms of n. 



could get a further constraint for the possible values of 
n using Eqn.JHJ). This equation will be taken instead as 
giving the current value of the snap as a function of rt 
(see fig.(|2Jl)|22- We have also plotted in the figure the 
error associated with so (dashed line), which grows as n 
for large n. The plot shows that with the current error of 
the kinematical parameters the method outlined here is 
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FIG. 2: Plot of the snap for f(R) given in Eqn.0 (solid 
curve) and the associated error (dashed curve) in terms of n. 



helpful in determining sq for the theory given by Eqn.J7|) 
only for n = 1. 

Let us next analyze an example of a theory involving 
two dimension- full parameters, given by |18| 
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f(R) =R + a\n—, 
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where // < 0. In this case it follows from Eqn.© and 
a < that 
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Figure © shows the permitted values for /j,/Rq in the 
case a < 0, which are between the horizontal axis and 
the solid curve, as well as the associated error j^. The 
plot shows that /i/Rq tends to a constant value (« 2.1) 
for large j3. Hence, the possible values of \i are restricted 
to \y\ < 1.2 x 1(T 41 uT 2 . 

In the case of the theory given in Ecm. HlUfl . the snap 
would be a function of a, /i, and of the remaining kine- 
matical parameters. 



Discussion 

A new method to restrict gravitational theories de- 
scribed by functions of the Ricci scalar has been intro- 



duced. It is based essentially in the assumption that 
normal matter composes the universe, the acceleration 
being caused by new gravitational dynamics in the low 
curvature regime, described by f{R). By imposing that 
the matter satisfy the weak energy condition, we obtain 
an inequality that constrain the parameters in the theory. 
We have shown by way of two examples how the method 
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FIG. 3: Plot of n/Ro (solid curve) and the associated error 
(dashed curve) in terms of a/Ro. The allowed values for fi 
are those below the curve. 



can be used, and how it conduces to restrictive limits on 
the parameters, having taken the error into account. We 
also obtained an equation that depends on the snap, the 
fourth derivative of the scale factor. Had we any mea- 
surements of so, this equation would furnish yet another 
condition on the parameters of the theory. Since the 
current value of the snap has not been determined yet, 
we take this equation as forecasting, for a given f(R), 
the current value of the snap. The method presented 
here could be combined with other approaches (such as 
avoidance of super- luminal propagation speed | l9| , com- 
patibility with the PPN limit [2£| , or those mentioned in 
the introduction) to restrict the f(R) theories that are 
being presented as candidates to model the acceleration 
of the universe. 
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